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ABSTRACT 

We present high and intermediate resolution spectroscopy of the X-ray nova Cen- 
taurus X-4 during its quiescent phase. Our analysis of the absorption features supports 
a K3 — K5 V spectral classification for the companion star, which contributes ~ 75 per 
cent of the total flux at Ha. Using the high resolution spectra we have measured the 
secondary star's rotational broadening to be usini = 43 ± 6 km s~^ and determined 
a binary mass ratio of g = 0.17 ± 0.06. Combining our results for K2 and q with the 
published limits for the binary inclination, we constrain the mass of the compact ob- 
ject and the secondary star to the ranges 0.49 < Mi < 2.49 M0 and 0.04 < M2 < 0.58 
Mq. a Doppler image of the Ha line shows emission coming from the secondary star, 
but no hotspot is present. We discuss the possible origins of this emission. 

Key words: accretion, accretion discs - binaries: close - stars: individual: Centaurus 
X-4 - X-rays: stars 



1 INTRODUCTION 

X-ray novae (XRNe) are a subclass of low mass X-ray bina- 
ries characterized by the exhibition of occasional outbursts 
with a recurrence time-scale of a few decades (see the reviews 
by van Paradijs & McClintock 1995; Tanaka & Shibazaki 
1996; van Paradijs 1998). During the outburst state the op- 
tical spectrum of the companion star is overwhelmed by the 
light of the X-ray irradiated accretion disc surrounding the 
compact object. During the long intervals of quiescence, the 
X-ray luminosity falls by at least 5 orders of magnitude, ren- 
dering the spectrum of the secondary star detectable. This 
provides an ideal opportunity to obtain the system parame- 
ters, constrain the physical properties of both the compact 
object and companion star and hence cast light on the na- 
ture and origin of these binary systems (see e.g. Charles 
1998). In this way it has been shown that ~ 75 per cent of 
the XRNe harbour black holes. 

Centaurus X-4 is a well-studied XRN. It was discovered 
during an X-ray outburst in 1969 (Conner, Evans & Belian 
1969), and a second outburst ten years later allowed the 
identification of its optical counterpart which had bright- 
ened by > 6 magnitudes (Canizares, McClintock & Grid- 
lady 1980). A Type I X-ray burst was detected during the 
decline phase of the outburst (Matsuoka et al. 1980), indi- 



cating that the primary in this case is an accreting neutron 
star. Photometric and spectroscopic studies of the optical 
counterpart in its quiescent state led to the determination 
of a 15.1 hr orbital period (Chevalier et al. 1989; McClintock 
& Remillard 1990) and a ~ 0.2 M© mass function (Cowley 
et al. 1988). The secondary star was identified as a K5 — K7 
evolved main sequence star. 

The Cen X-4 quiescent light curves show a double- 
humped modulation with a full amplitude of ~ 0.15 mag 
(Cowley et al. 1988; Chevalier et al. 1989; McChntock & 
Remillard 1990) which is associated with the tidal distor- 
tion of the Roche-lobe-filling secondary star. McClintock & 
Remillard (1990), based on the phasing of the H/3 emission 
line radial velocity curve, suggested the presence of an ac- 
cretion disc hotspot to explain the unequal maxima and un- 
usual different depths of the minima in the V and / light 
curves. Star spots or X-ray heating of the K-star were also 
invoked to explain the long term changes observed in the 
light curves by Chevalier et al. (1989). Shahbaz, Naylor & 
Charles (1993) observed Cen X-4 in the IR, where the con- 
tamination of the ellipsoidal light curve by the accretion 
disc is minimized. Their fit of the _ff-band light curve with 
an ellipsoidal model led to the determination of a low incli- 
nation (31° — 54°) and a primary mass lying in the range 
0.5 - 2.1 Mq. 
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In this paper we refine the system parameters and 
ephemeris of Cen X-4. From our higher resolution spectra 
we are able to determine the rotational broadening of the 
companion star for the first time. Finally, we investigate the 
Ha Doppler map of the system during the quiescent state. 



2 OBSERVATIONS AND DATA REDUCTION 

Centaurus X-4 in its quiescent state was observed during 
three runs. The first was from 1993 June 24-27 when spec- 
tra were acquired using the RGO spectrograph attached to 
the 3.9-m Anglo- AustraUan Telescope (AAT). A 1200 fine 
mm~^ grating centred on Ha and a 1.2 arcsec slit width 
yielded a dispersion of 0.8 A pix~^ and a spectral resolution 
of 74 km s"^ (FWHM). On 1994 March 6-8 we employed the 
EMMI spectrograph on ESO's New Technology Telescope 
(NTT) on La Silla. Use of grating # 6 and a 1.2 arcsec slit 
yielded a dispersion of 0.31 A pix~^ and a spectral resolu- 
tion of 45 km s~^ . Our final run was with the 4.2-m William 
Herschel Telescope (WHT) on La Palma. Spectra were ac- 
quired during the nights of 1994 July 9/10 and 13/14 using 
the red arm of ISIS, a 1200R grating and a 0.8 arcsec slit 
which provided a spectral resolution of 25 km s~^. In addi- 
tion to the spectra of Cen X-4, on the night of July 9 radial 
velocity standard stars of spectral types K3 HI (HR 5536) , 
K4 HI (HR 6159) and K5 HI (HR 5622) were acquired using 
the same setting as for the target. Table |^ provides a more 
detailed journal of the Cen X-4 observations. 

All the runs were photometric, with seeing less than 
1.5 arcsec at all times, and < 1 arcsec occasionally. Copper- 
Neon arc lamp spectra were taken every hour. The spectra 
were extracted in a standard way using irafQ. The pixel-to- 
wavelength calibration was derived from third-order polyno- 
mial fits to 12 arc lines, giving a rms scatter < 0.03 A. The 
calibration curves were interpolated in time and the accu- 
racy of the whole process was checked by cross-correlating 
the sky spectra, which was found to be within 10 km s^^. 
The spectra were re-binned into a constant velocity scale 
and normalized by dividing with the result of fitting a low 
order spline to the continuum after masking the emission 
lines. 



3 RADIAL VELOCITY CURVE AND ORBITAL 
EPHEMERIS 

Radial velocities were extracted by cross-correlating each 
target spectrum with the spectra of the three template stars 
in the range AA6387 — 6513. We performed least-squares 
sine fits to our radial velocity data, fixing the orbital pe- 
riod at P = 0.629063 d (McCfintock & Remillard 1990). 
The resultant dynamical parameters are summarized in Ta- 
ble 1^. We adopt the K4 HI parameters of the radial ve- 
locity curve because it provided the lowest reduced x^: K2 
= 150.0 ± 1.1 km s-\ 72 = 184.0 ± 1.9 km s"^ and 
To = HJD 2449163.934 ± 0.001. This zero phase is defined 



* IRAF is distributed by the National Optical Astronomy 
Observatories. 



as inferior conjunction of the secondary star. The combina- 
tion of our To determination with the value reported by Mc- 
Clintock & Remillard (1990) enables us to refine the orbital 
period to P = 0.6290496 ± 0.0000021 d. Fig. displays the 
radial velocity curve. Although our velocity semi-amplitude 
K2 agrees with the values presented by Cowley et al. (1988) 
and McClintock & Remillard (1990), we find a systemic ve- 
locity intermediate to these reported earlier This discrep- 
ancy could be due to the use of an erroneous heliocentric 
radial velocity for the template stars used by these authors. 
Unfortunately we have not found any updating of them. 



4 SPECTRAL CLASSIFICATION AND 
ROTATIONAL BROADENING 

We determined the spectral type and rotational broadening 
of the secondary by subtracting a set of templates from the 
Doppler-corrected average spectrum of the target until the 
lowest residual was obtained (see e.g. Marsh, Robinson & 
Wood 1994). 

Regarding the spectral classification we have con- 
strained this analysis to these spectra acquired with the 
AAT because of their higher signal-to-noise ratio. The target 
spectra were summed in the rest frame of the template star 
HR 6159 by using the radial velocity parameters of Sect. 3. 
Different weights were assigned to different target spectra in 
order to maximize the signal-to-noise of the sum. In addi- 
tion to the K Ill-type standards acquired with the WHT, we 
have employed a set of K V-type standards observed with the 
Isaac Newton Telescope (INT) in 1993 October (see Casares 
et al. 1996). These spectra were Doppler-shifted to the rest 
frame of HR 6159 and degraded appropriately to match the 
broadening of our AAT observations. This value was found 
by subtracting difi^erent broadened copies of the WHT and 
INT spectral templates from the Doppler-corrected average 
of Cen X-4 and performing a test on the residuals. 

The templates were multiplied by a factor < / < 1, 
representing the fractional contribution of light from the 
secondary star, and subtracted from the target average. 
Then a test on the residuals was performed in the range 
AA6387 — 6513. The optimal values of / are provided by 
minimizing , with quoted uncertainties corresponding to 
Xmin+1 (Lampton, Margon & Bowyer 1976). The minimiza- 
tion of xi does not give a single minimum (see Table ^) and 
the secondary star in Cen X-4 is most probably of spectral 
type K3 — K5 contributing ~ 75 per cent to the total flux of 
the system at Hq. A similar result was obtained when ap- 
plying this technique to the NTT and WHT data subsets. 
Because HD 29697 is an active star (Casares et al. 1996), this 
might hint that the companion star in Cen X-4 is also ac- 
tive. Fig. 1^ depicts the normalized Doppler-shifted average of 
Cen X-4 before and after subtraction of the K5 V template. 
The Cen X-4 spectrum shows Ha and He I A6678 emission 
lines. Li I A6708 in absorption becomes prominent after the 
subtraction. The detection of a strong Li I A6708 line in the 
late-type companion is an unexpected result since the initial 

t These values are: K2 = 152.0 ± 10 km s"!, 7 = 234 ± 8 km 
s-i (Cowley et al. 1988). K2 = 146 ± 12 km s"!, 7 = 137 ± 17 
km s-i (McClintock & Remillard 1990). 
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Li I content should be quickly depleted by means of convec- 
tive mixing in the star and mass transfer to the compact 
object. The reader should refer to Martin et al. (1994) for a 
detailed discussion of the significance of the anomalous high 
abundance of Li I in Cen X-4 and other XRNe. 

For the rotational broadening calculation we have se- 
lected a subset of our data which has the higher spectral 
resolution. This consists of the target and template spectra 
acquired with the WHT. The spectra were Doppler-shifted 
to the rest frame of the template HR 6159 and then the tar- 
get spectra were summed with different weights to optimize 
the final signal-to-noise ratio. In order to correctly estimate 
vsini we simulate in our template spectra the drift due to 
the orbital motion of the secondary during the exposures: 
we created 8 versions of each WHT template which were 
smeared in the same way as the target spectra by convo- 
lution with a rectangular profile of the appropriate width. 
This is given by: 

O-smcar = * cxp "p -^^"2 | COS {2TVip) \ 

where tcxp is the length of the exposure and tp the orbital 
phase. 

The smeared template copies were then summed us- 
ing identical weights as for the Cen X-4 spectra. Subse- 
quently, we broadened the three summed templates from 
1 to 100 km s~^, in steps of 1 km s~^, and subtracted 
the resultant broadened spectra (as above) from the tar- 
get Doppler-corrected average. The summed templates were 
broadened through convolution with the rotational profile of 
Gray (1992). We used linearized limb-darkening coefficients 
of fi = 0.69, n = 0.73 and /i = 0.75 appropriate for a K3 
III, K5 III and K7 III star, respectively (Al-Naimiy 1978; 
Wade & Rucinski 1985). The results from the xt minimiza- 
tion are listed in Table W, with quoted uncertainties on v sin i 
corresponding to Xmin+l- Because the minimum xi for the 
output residuals is acceptable (~ 1) for the three templates, 
we choose for the rotational broadening the mean of the re- 
sults obtained from the analysis: vsini = 43 ±6 km s"\ To 
check the systematic errors introduced by the choice of /i 
(which is only suitable for the continuum; Collins & Truax 
1995), we have allowed fi to vary in the range 0-1. This leads 
to 5 per cent (~ 3 km s~^) changes in the resulting value 
of vsini, and therefore it is the statistical noise and not the 
uncertainty in fi which limits our accuracy. Unfortunately, 
we cannot make use of the higher signal-to-noise spectra 
acquired with the AAT and NTT to refine the rotational 
broadening. The absorption features of the former are artifi- 
cially broadened by the instrumental resolution (74 km s~^), 
whereas the latter are dominated by the instrumental reso- 
lution (45 km s~^) and the orbital smearing. In particular, 
the total broadening (calculated by adding in quadrature 
the instrumental resolution and (Tsmear) ranges from 25 to 
37 km s"^ for the WHT spectra and 47 to 61 km s"^ for the 
NTT spectra. 



5 THE BINARY SYSTEM PARAMETERS 

Assuming that the companion star is synchronized with the 
binary motion and fills its Roche lobe we can calculate the 
meiss ratio {q = M2/M1) through the expression (see e.g. 
Wade & Home 1988): 



vsini = 0.462 K2 q^^^(l + qf'^ 

From the values of v sin i and K2, we derived q — 0.17 ±0.06, 
which implies a primary velocity semi-amplitude of Ki = 
qK2 = 26±9 km s~^. On the other hand, our values of K2 — 
150.0±1.1 km s~^ and P = 0.6290496±0. 0000021 d imply a 
mass function of /(M) = 0.220 ± 0.005 Mq. Combining the 
values for f{M) and q, the masses of the compact object and 
the companion star are: Mi = (0.30 ± 0.04) sin~^ i Mq and 
M2 = (0.05±0.03) sin"^ z M©. We used the constraint upon 
the inclination from the IR observations of Shahbaz et al. 
(1993). By fitting the _ff-band light curve with an ellipsoidal 
model these authors found the system inclination to be in 
the range 31° < i < 54° (90 per cent confidence). Taking 
the extreme limits for i we obtain ranges for the masses of: 
0.49 < Ml < 2.49 Mq and 0.04 < M2 < 0.58 Mq. The 
companion star of Cen X-4 is of spectral type K3 — K5 V. 
Main-sequence stars of these spectral types have masses in 
the range 0.74 — 0.68 Mq. The disagreement with the masses 
obtained above is easily explained if the secondary star is 
undermassive (i.e. evolved) for its spectral type as proposed 
by Shahbaz et al. (1993). 



6 DOPPLER TOMOGRAPHY 

The Ha emission line has a FWHM of 640 km s"'^ and a 
mean equivalent width (EW) of 35 ± 7 A, with the er- 
ror accounting for the intrinsic line variability during the 
three runs. The line profile has a double-peaked structure, 
the classic signature of an accretion disc, and shows peri- 
odic changes in the height of the two peaks by the action 
of an S-wave emission component. We have made use of 
the Doppler tomography technique (Marsh & Horne 1988) 
on the orbitally resolved line profiles to reconstruct the Ha 
brightness distribution and identify the location of the S- 
wave in the binary system. In order to construct the Doppler 
map, the 25 spectra were re-binned into 82 pixels of constant 
velocity scale (36.6 km s~^) and the continua removed. The 
derived Doppler image computed with the maximum en- 
tropy method (Marsh & Horne 1988) is shown in the bottom 
panel of Fig. H. This image should be interpreted as a one 
year time-averaged map since the data were acquired over 
a year. We have overplotted the theoretical path of the gas 
stream and the Roche lobe of the secondary star and the gas 
stream trajectory for K2 — 150 km s~^ and q — 0.17. The 
Doppler image shows a region of enhanced intensity at the 
position of the secondary star, superposed on a weak ring- 
like structure, which represents the emission arising from a 
rotating accretion disc. There is no evidence for emission 
from the gas stream or the hotspot (the stream/disc impact 
region) . We think unlikely that this absence is an artifact of 
using spectra spread over a year since a very similar result 
is obtained when using only the AAT data subset which ex- 
pands over 5 consecutive orbital cycles. Individual Doppler 
maps for the WHT and NTT data subsets could not be com- 
puted because of scarce phase coverage. However, we note 
that the S-wave component from the secondary is clearly 
detected in all the individual line profiles. 
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7 DISCUSSION 

In order to constrain the possible causes of the Ha emission 
arising from the secondary, we have extracted the Ha S-wave 
component from the line profiles as described by Casares et 
al. (1997). The procedure applied is illustrated in Fig. ^ 
First, we simulate the accretion disc by averaging all the 
spectra in the laboratory rest frame. Because the resultant 
profile showed the blue peak significantly stronger than the 
red one (clearly an artifact of the S-wave emission due to the 
gaps in our phase coverage), we built a symmetric double- 
peaked profile by folding a copy of the red part of the emis- 
sion line profile (spectrum a in Fig. Next, twenty five 
versions of this 'accretion disc' profile were shifted to the 
rest frame of the secondary star, at the times of our indi- 
vidual spectra, and then averaged using the same weights 
(spectrum h in Fig. ^). Finally, the resultant profile was re- 
moved from the Doppler-corrected average of Cen X-4 with 
the K5 V template spectrum subtracted, to be left with the 
narrow Ha component with EW(Ha) = 1.7±0.5 A (spectrum 
d in Fig. In an attempt to find secular variations in the 
strength of the Ha S-wave component, we extracted it (as 
above) from each individual data subset. It was found that 
the values of EW(Ha) were comparable within the errors to 
the one obtained using the complete baseline. 

To compare EW(Ha) with observations of rapidly ro- 
tating stars, we convert this value into Ha line surface flux 
at the star (-Fka). Fna is derived using the relationship given 
by Soderblom et al. (1993): 

log Fhc = log EW(Ha) + 0.113(5 - V)l 
-1.188(B - V)o + 7.487 

Using [B - y)(,=-f 1.25 mag (Chevalier et al. 1989) and an 
equivalent width (corrected by the disc's contribution) of 
EW(Ha)= 1.7/0.75 = 2.3 A, we obtain Fhq = 3.5 x lO'' erg 
cm~^ s~^. The dependence of the chromospheric activity 
on rotation is given by the relation between the fraction 
of the stellar luminosity emitted in the Ha line (7?hq = 
FYia/o-T^s) and the Rossby number _Ro = P/tc (Noyes et al. 
1984), where P is the orbital period of the star and Tc is the 
convective turnover time. Adopting T^s = 3600 K (Shahbaz 
et al. 1993) and Tc ~ 25 d, as calculated using eq. 4 of 
Noyes et al. (1984), we derive log_RHa = —3.4 and logi?o = 
— 1.6. These values can be compared with observations of 
rapidly rotating (P > 0.25 d) chromospherically active stars 
in the Pleiades (Soderblom et al. 1993; fig. 19a). According 
to fig. 19a, saturation of chromospheric activity would occur 
for P < 2d or wsini > 15 km s~^. Indeed the companion of 
Cen X-4, with P = 0.63 d and vsini — 43 km s~^, lies in 
the region of activity saturation (logPna — —3.7), but with 
Rua slightly higher than expected. 

We now examine the possibility that the Ha compo- 
nent is caused by the incidence of X-rays on the inner face 
of the secondary. The observations of Cen X-4 have shown 
that its quiescent X-ray luminosity is variable on time-scales 
from days (Campana et al. 1997) to years (van Paradijs et 
al. 1987; Rutledge et al. 2001). We adopt the 1994 intrinsic 
X-ray luminosity of 2.4 x 10^^ erg s"^ (0.5 - 10 keV) ob- 
tained with ASCA (Asai et al. 1996). Using this value and 
neglecting any shielding by the accretion disc, the X-ray flux 
irradiation of the secondary star would be Fx = 3 x 10* erg 
cm~^ s~^, where we have used a binary separation of 3.6 R© 



(assuming Mi = 1.4 M©). At these flux levels, the secondary 
would be receiving sufficient energy to power the narrow Ha 
component. Therefore we can only conclude that the Ha 
emission arising from the secondary star may be either due 
to chromospheric activity, and/or powered by X-ray heating 
on the secondary. 

Ha emission with origin in the secondary star has also 
been found in the quiescent systems Nova Muscae 1991 
(P = 10.4 hr, vsini = 106 km s'^ Casares et al. 1997) 
and probably in A0620-00 (P = 7.8 hr, vsmi^ 83 km s"^ 
Marsh et al. 1994), V404 Cygni (P = 155.3 hr, vsmi = 39 
km s~^; Casares & Charles 1994; Casares 1996), Nova Ophi- 
uchi 1977 (P = 12.5 hr, -usini < 79 km s"^ Harlaftis et al. 
1997) and Nova Scorpu 1994 (P = 62.9 hr, -usini = 86 km 
s~^; Shahbaz et al. 1999). The high rotational broadening 
measured suggests that coronal activity may exist in these 
XRNe. Unfortunately, the poor phase coverage in the case 
of Nova Sco 1994 and Nova Oph 1977, the lack of EW val- 
ues of the Ha component in A0620-00 and V404 Cyg, and 
the fact that we only have an upper limit to the X-ray flux 
of Nova Mus 1991, prevent us from doing a more detailed 
analysis of the origin of the Ha emission in these systems. 

The contribution to the optical flux due to a hotspot 
is required in order to explain and model the optical and 
IR light curves of Cen X-4 (McChntock & Remillard 1990; 
Shahbaz et al. 1993). The presence of such hotspot was 
strongly supported by the phasing of the 1987 H/3 emission 
line radial velocity curve (McClintock & Remillard 1990). 
Therefore its absence in our Ha tomogram is a puzzling re- 
sult. A possible explanation would be that the hotspot is 
fainter or comparable in brightness to the accretion disc at 
Ha and brighter at shorter wavelengths (higher members of 
the Balmer series or high excitation lines) . A second possibil- 
ity is that the secondary star in Cen X-4 might under-fiU its 
Roche lobe by a few percent of its radius and mass transfer 
is unstable. 

The acquisition of new orbitally resolved high-resolution 
data will bring substantial progress in the knowledge of Cen 
X-4. A simultaneous and uniform orbital phase coverage of 
the binary orbit may reveal whether or not the origin of the 
secondary star's Ha emission is only chromospheric or due to 
X-ray heating. The absence of hotspot at Ha in XRNe in qui- 
escence has also been seen in Nova Persei 1992 (P = 5.09 hr; 
Harlaftis et al. 1999), Nova Mus 1991 (Casares et al. 1997) 
and maybe in Nova Oph 1977 (Harlaftis et al. 1997). This 
result suggests that it may be a common phenomena among 
X-ray transients. Unfortunately very little data following 
an XRN during its quiescent state are presently available. 
Therefore observations at different epochs of quiescence are 
required to show whether this is a transient phenomenon or 
a temperature effect of the hotspot. 



8 SUMMARY 

We have presented spectroscopic observations of the XRN 
Cen X-4 during quiescence. A test applied to the resid- 
uals obtained by subtracting different template stars from 
the Cen X-4 spectrum supports a secondary star of spectral 
type K3 — K5 V contributing ~ 75 per cent of the observed 
flux at Ha. Using the high-resolution spectra, we have deter- 
mined the rotational broadening of the companion star to be 
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V sin i = 43±6 km . This has provided a rehable measure- 
ment of the mass ratio in Cen X-4: q — 0.17±0.06. We obtain 
a mass function of 0.220 ± 0.005 Mq which, when combined 
with the mass ratio and binary incMnation constraints, gives 
a range of 0.49 < Mi < 2.49 Mq and 0.04 < Ah < 0.58 Mq 
for the neutron star and secondary star masses. This result 
reinforces the possibiUty that the secondary is evolved for its 
spectral type. Furthermore, we have revised the ephemeris 
of Cen X-4 for the inferior conjunction of the secondary 
to To = HJD 2449163.934(1) + 0.6290496(21) x N d. Our 
Doppler tomogram of the Ha line is dominated by a strong 
emission arising from the companion star but no hotspot 
contribution is observed. The Ha flux derived from the S- 
wave component is consistent with the value expected for a 
chromospherically active and/or X-ray heated secondary. 
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Table 1. Journal of Observations. 



Date Telescope No. spectra Exp. time A range Dispersion 

(s) (A) (A pix-i) 

24-27 Jun 1993 AAT 11 1800 6155-6968 0.80 

6 Mar 1994 NTT 3 2400 6175-6815 0.31 

7 Mar 1994 " 2 1800-1900 

8 Mar 1994 " 1 3300 

9 Jul 1994 WHT 2 1800 6360-6767 0.40 

10 Jul 1994 " 2 1200-1800 

13 Jul 1994 " 2 1800 

14 Jul 1994 " 2 " 6307-6714 



Table 2. Radial velocity parameters (errors at 1<t). 



Template Spectral 7 7^ Tq xl 

Type (km s"!) (km s"!) (km s"!) (4-2249163) (d.o.f=22) 

HR 5536 K3 111 207.4±0.9 187.0 ± 2.9 152.2 ± 1.0 0.934 ± IxlO^^ 7.1 
HR6159 K4 III 180.9±1.0 184.0 ± 1.9 150.0 ± 1.1 0.934 ± 1x10"^ 4.5 
HR5622 K5 III 200.5±0.8 185.4 ± 1.7 151.0 ± 1.0 0.934 ± IxlQ-^ 6.8 

" 72 : Systemic velocity after adding the radial velocity of the standard star found in SIMBAD. 



Table 3. Spectral classification. 



Template 


Spectral 




/ 




Type 


(d.o.f.=159) 




HR 5536 


K3 HI 


1.92 


0.72 ± 0.03 


HR 6159 


K4 HI 


1.97 


0.76 ± 0.03 


HR 5622 


K5 HI 


2.11 


0.64 ± 0.02 


HD 184467 


K2 V 


1.75 


1.00 ± 0.03 


HD 29697 


K3 V 


1.30 


0.72 ± 0.02 


HD 154712 A 


K4 V 


1.48 


0.76 it 0.03 


61 Cyg A 


K5 V 


1.34 


0.75 ± 0.03 


61 Cyg B 


K7 V 


1.58 


0.69 ± 0.02 



Table 4. Rotational broadening. 



2 

mi 

Type (km s"!) (d.o.f.=321) 



Template Spectral n v sin i X'min v 



HR 5622 K3 HI 0.69 42.5 ± 6.0 0.97 
HR 6159 K4 111 0.73 44.2 ± 5.4 0.97 
HR 5536 K5 HI 0.77 42.0 ± 5.8 1.00 
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Figure 1. Radial velocities of the secondary star in Cen X-4 
folded on the epheineris of Sect. 3.1. Observations were performed 
in 1993 June (solid squares), 1994 March (squares) and 1994 July 
(diamonds). The best sine- wave fit is also shown. 
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Figure 2. Optimal subtraction of a K5 V template star (bottom 
spectrum) from the Doppler-corrected average of Cen X-4 in the 
rest frame of the secondary (middle spectrum) . The residual of the 
subtraction (top spectrum) represents the relative contribution of 
the accretion disc. Arbitrary vertical offsets have been applied to 
the spectra for the sake of clarity. 
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Figure 3. Doppler image of the Ha emission (bottom panel). The 
data, represented in the form of a trailed spectra, are depicted in 
the top panel. Empty strips represent gaps in the phase coverage. 
For the sake of clarity, the same cycle has been plotted twice. The 
Roche lobe of the secondary star, the predicted velocities of the 
gas stream (lower curve) and of the disc along the stream (upper 
curve) arc plotted for K2 = 150 km and q = 0.17. Distances 
in multiples of O.IRli arc marked along both curves with open 
circles. The centre of mass of the system is denoted by a cross. 
Note that the Doppler image has been computed from the origi- 
nal data without subtracting the spectrum of the secondary star, 
improper subtraction of which can lead to an erroneous indication 
of emission from the secondary (see e.g. Marsh et al. 1994). 
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Figure 4. Dctrcnding the narrow Hq coiiiponciit. The simulated 
symmetric Hq profile (a) has been smeared according to the or- 
bital phases of the individual spectra (b). Removal of (b) from 
the Doppler-corrected average spectrum of Cen X-4 in the rest 
frame of the secondciry after subtraction of the K5 V template 
(c) yields the narrow Ha component (d). Arbitrary vertical off- 
sets have been applied to the spectra for the sake of clarity. 
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